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ABSTRACT 


This thesis develops a block diagram representation of a salient 
pole synchronous generator. This representation in its final form 
provides an accurate means for analyzing the machine by analog com- 
puter methods. 


Linear volt-ampere equations describing the synchronous 
alternator are developed from the basic equations for the various 
machine circuits by means of matrix algebra. Using these linear 
volt-ampere equations and making certain simplifying assumptions, 

a single time constant computer block diagram which represents the 
alternator and its connected load is developed. The nonlinear effects 
of saturation are introduced at a single point in this block diagram. 


The validity of the representation is shown by a comparison of 
steady state and transient phenomena as obtained experimentally from 
the alternator and analytically from the computer representation of 
the alternator. The steady state results show exact'agreement of the 
open circuit condition and only slight divergence in the upper regions 
of saturation for the experimental and computer curves describing 
operation at both zero and unity power factors with rated armature 
current. 


Transient phenomena are measured by application of both load 
steps and field excitation steps for three widely separated conditions 
of loading. Agreement between experimental and computer transients 
is excellent for the load changes and acceptable for the field excita- 
tion changes. 
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Conclusions 


(1) 
(2) 
(3) 


The block diagram developed is a valid representation of the 
synchronous alternator. 


The final representation is ina form readily adaptible to computer 
presentation. 


The experimental and computer steady state characteristics curves 
diverge somewhat at high values of field excitation for both unity 
and zero power factor loads. This divergence will decrease for 
intermediate values of loading. 


Recommendations 


The recommendation is made that the following points should be 


investigated further: 


(1) 
(2) 


Methods to reduce the steady state errors at high values of field 
excitation with both zero and unity power factor loads. 


Plausibility of including the effects of armature transients and 
amortisseur subtransients in the alternator block diagram. 


(3) A method to properly account for speed variations in the repre- 


(4) 


sentation. 


Study of inclusion of unbalanced loads in the generator-load repre- 
sentation. 
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INTRODUCTION 


The transient performance of synchronous machines has been 
the subject of many careful investigations which have resulted in a. 
considerable amount of data representing operation under all types of 
conditions. Most of these investigations have been based upon linear 
theory in order to keep the analysis as clear and uncomplicated as 
possible. When included in these analyses, nonlinear effects were 
introduced by adjustment of machine parameters at a particular con- 


dition of operation and were valid only at that condition. 


A more recent investigation (11)(12) by the Energy Conversion 
Group, Massachusetts Institute of Technology, adopted a somewhat 
different approach. This investigation was primarily concerned with 
the voltage regulatory problem of an alternating current generator 
system. Since the various dynamic problems associated with regula- 
tion, such as load changes, faults, and switching, as well as the many 
components of the system itself, are nonlinear, a completely linear 
analysis was impossible. Consequently, one of the purposes of this 
investigation was the development of a block diagram representation 
of the voltage regulatory system, including nonlinear effects, which 


would be suitable for a study of the problem using an analogue computer. 


In the portion of the study concerned with the synchronous 
alternator the Energy Conversion Group was successful in obtaining 
a computer representation of the alternator which proved to be valid 


with a zero power factor load. 








The purpose of this report is to present the steps which were 
followed in the development of a dynamic computer block diagram 
representation of a synchronous alternator, including nonlinear char- 
acteristics, which is valid for any type of transient disturbance and 


for loads of power factors from zero to one inclusive. 





II 


PROCEDURE 


The procedure followed in developing the computer block diagram 
representation of a synchronous alternator and testing its validity can 


be divided into the following steps: 


Development of the Linear Volt-Ampere Equations Describing the 
Synchronous Alternator 

In this portion of the report the electrical equations of motion for 
the stator and rotor were written and, by a series of transformations, 
were changed to direct and quadrature two-axis components (10). The 
equations for the synchronous machine were first developed in motor 
form and were then changed to the generator form. Throughout this 
development matrix algebra was used to provide a compactness of form 
with the transformations following a path from three-phase components 
to Lyon's symmetrical components, to forward and backward rotating 
components, and finally to Park's direct and quadrature two-axis com- 
ponents. In order to obtain a final set of volt-ampere relationships 
which are linear differential equations with constant coefficients, it 


was necessary to assume constant speed operation. 


Appendix A is a more thorough description of this portion of the 


report. 


Development of a Computer Block Diagram Representation of the 


Synchronous Alternator and its Connected Load 


A block diagram was developed from the volt-ampere equations 





describing the synchronous alternator. In order to include a load 
representation in this block diagram a vector diagram describing the 
machine in two-axis components was constructed. Since the block 
diagram developed was not suitable for direct application to an analog 
computer, it was necessary to change this block diagram to computer 
form and include a valid representation of the nonlinear characteristics 


of the synchronous alternator. 


Appendix B describes the development of a computer represen- 


tation more completely. 


Experimental Determination of Machine Parameters and Transients 
Phenomena 

After choosing a suitable synchronous generator, the machine 
parameters were determined from the following tests and measure- 
ments: slip test, open circuit test, short circuit test, and a measure- 
ment of the field time constant. This experimental work provided the 
necessary machine parameters which were then converted to per unit 
values in the reciprocal system described by reference 8. Resulting 
from this conversion were the per unit values necessary to complete 
the computer block diagram representation. These are listed below. 


R, = per unit field resistance including field rheostat 
and external circuit resistances, i 


X nd - per unit magnetizing reactance, 
Xe = per unit field leakage reactance, 


X,, = per unit armature leakage reactance, 


al 
and the nonlinear relation between resultant open circuit field 
current, ip , and "air gap" field current, is 

oc ag 





The transient characteristics of the synchronous generator were 
determined at three widely separated load conditions: pure inductive 
load (pF = 0), pure resistive load (pF = 1), and a combination inductive- 
resistive load (pF = 0.752). The choice of load also determined the 
remaining two blocks, Y K (9) < +1 and K , (8), in the load portion of 
the computer representation. For each of these load conditions the 
dynamic effects upon armature terminal voltage, armature current, 
and field voltage were determined by oscillograms when making the 
following setp changes: no load to full load, full load to no load, and 


a field voltage step. 


The final computer block diagram representation was based 
upon the assumption that the machine speed remained constant. During 
the determination of the transient phenomena, the machine speed 
varied upon application of the various step changes, the amount of 
variation depending upon the load condition. An oscillogram was taken 
of this speed transient in order to determine the effect of the speed 


variation upon the other transients measured. 


Appendix C is a more complete description of the experimental 


determination of machine parameters and transient phenomena. 


Computer Determination of Transient Phenomena 


The alternator representation obtained should be suitable for 
analytic study and solution by analog computer methods, and was 
developed with this object in mind. In order to test the validity of the 
alternator representation, the solution of the block diagram should be 
checked with transient and steady-state values obtained by direct ex- 


perimentation on the test machine. An initial computer test was, 





therefore, to obtain characteristic curves of steady-state relationships 
between field excitation and generator voltage under various conditions 
of loading. These curves should compare with the saturation curves 
experimentally obtained. The computer study of transient phenomena 
was made duplicating the field excitation and load conditions which 


were used in the experimental transient study. 


Appendix D describes the computer analysis more completely. 


III 


RESULTS 


Figure I presents the block diagram representation of the 
salient pole seis non oui generator. This representation was used 
to determine various steady state and transient terminal voltage 
responses from an analog computer study. An experimental study 
was also made duplicating the conditions used for obtaining the com- 
puter responses. This section of the thesis correlates the results 
of the experimental and computer work and proves the validity of 


the developed block diagram. 


Figure XI presents a comparison of the saturation curves 
obtained by both experimental and computer study at no-load and at 


zero and unity-power factors with rated armature current. 


Figures II through X present a comparison of the transient 
responses of the terminal voltage for changes in load and field exci- 
tation as obtained by both computer and experimental studies. The 
initial and final voltages in the case of each transient were adjusted 
on the experimental curves to coincide with the initial and final 
voltages obtained from the computer curves. This was made neces- 
sary because the computer study was carried out assuming constant 
speed, whereas actually the experimental transients were accom- 


panied by generator speed changes. 


The computer representation was developed by neglecting 


armature transients and subtransient effects due to the amortisseur 





windings. The computer transients reflect this oniission in the 
theoretical development by a sudden change at the point where the 
transient was initiated. These transient and subtransient reactance 
changes are present in the experimental curves near the point of 


transient initiation; however, they die out very rapidly. 


Table I is a comparative tabulation of the time constants 
obtained from the experimental and computer transients as measured 
when neglecting those transients due to armature reactance and 


amortisseur subtransient reactance. 





TABLE I 


TIME CONSTANT COMPARISON 


Experimental 

Power factor = 0 

Full load to no load 0.10 

No load to full load 0.125 

Field step 0.158 
Power factor = 1.0 

Full load to no load 0.125 

No load to full load 0.13 

Field step 0.148 
Power factor = 0.752 

Full load to no load 0.11 

No load to full load 0.14 

Field step 0.145 


Computer 


0.105 


0.145 


0.175 


0.145 


0.295 
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IV 


DISCUSSION OF RESULTS 


Steady State Relationships 


The no load saturation curves obtained by computer and 
experimental methods compared exactly over the entire range of field 
excitation used. This is as might be expected since the circuit at no 
load in the block diagram representation is completely defined by fixed 
quantities, i.e., the nonlinear characteristic, the field voltage, and 
field impedance. Thus, the no load correlation merely verifies the 


construction of the nonlinear characteristic curve. 


The zero power factor saturation curves for both the computer 
and experimental studies were found to compare exactly as field cur- 
rent was increased until a point was reached slightly above rated ar- 
mature voltage. At that point the computer curve began falling below 
the experimental curve with slowly increasing divergence. The analy- 
tical (computer) representation was developed assuming that the arma- 
ture leakage reactance was a constant quantity equal to the Potier 
reactance (see Appendix C). Thus, the same Potier triangle moved 
parallel to itself would fit between the open circuit and zero power 


factor curves throughout the range of field excitation. 


The zero power factor saturation curve does not have exactly 
the same shape as the open circuit characteristic shifted downward 
and to the right by constant amounts. The principal reason for this 


discrepancy is that the field leakage flux is not a constant under varying 


15 





conditions of loading. At zero power factor and with higher field 
excitations than those required for rated armature voltage, the field 
current is greater than for the same armature voltage on open cir- 
cuit. This causes a higher field leakage flux and increased saturation 
due to the main air gap flux and the field leakage flux. This effectively 
causes the armature leakage reactance drop to be somewhat less than 
the Potier reactance. In general, the presence of field leakage ina 
salient pole machine causes the open circuit and load characteristic 
curves to draw near each other vertically for high values of saturation. 
This subject is covered with considerable clarity by Fitzgerald and 


Kingsley (see reference 4). 


The correlation between the unity power factor saturation curves 
obtained experimentally and by computer is fairly close, agreeing 
exactly in the lower linear regions and diverging somewhat with in- 
creasing values of saturation. The computer curve gives higher arma- 
ture voltage than the experimental curve at the same field excitation 
when in the saturated region. This discrepancy is largely due to the 
assumption of no saturation of magnetic material in the quadrature 
axis. This assumption is not entirely correct, but because the inter- 
polar air gap spaces have such a high proportion of air in the magnetic 
circuit as compared to that in the direct axis magnetic circuit, the 
value of quadrature axis synchronous reactance, Xy was assumed 
constant and unsaturated. Actually, some degree of saturation does 
exist at higher values of field excitation due to the presence of iron in 
the quadrature axis magnetic circuit, although not as much as in the 
direct axis. The effect of quadrature axis saturation is to reduce the 


value of X_ from its unsaturated value. Since the direct axis terminal 
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voltage, ey is directly proportional to Ka’ saturation has the effect 
of reducing ea and, consequently, the terminal voltage. Thus, the 
experimental saturation curve for unity power factor would rightly be 


expected to be lower than the corresponding computer curve. 


It should be noted that while the computer characteristic for 
unity power factor is higher than the unity power factor experimental 
curve, the reverse is true for zero power factor loads as earlier ex- 
plained. This is evident from a study of figure XI of the results. It 
therefore appears quite likely, that at some intermediate power factor 
between zero and unity, the variation of the analytic steady state values 
from those experimentally determined would be zero. For sucha 
power factor the effect of decreasing field reactance with saturation 


would be cancelled by the decreasing value of Xa with saturation. 


It would be possible to reduce the discrepancy between computer 
and experimental curves at unity power factor by inserting a satura- 
tion factor at * in the block diagram of figure I of the section on 
results to account for the discrepancy which presently exists. Further 
compensation to account for the small discrepancy between the com- 
puter and experimental zero power factor curves could also be accom- 
plished by insertion of a nonlinear factor at the location of X The 
resulting block diagram would give more accurate analytic results but 
would involve three nonlinear elements hic are characteristic of a 
single machine. The relative simplicity of the representation developed 
in this thesis would have to be sacrificed to some degree in order to 
obtain a closer correlation between the experimental and analytic 


results. 
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Transient Relationships 


Figures 11 through X of the section on results graphically compares 
the transient phenomena observed for the experimental and computer 
studies. The comparisons for full load to no load and for no load to 
full load are quite excellent. The transients in the terminal voltage 
due to a field step are also good, although not agreeing as closely as 


for the load changes. 


The experimental transient is really a three time constant 
curve, one due to field transients, another due to armature transients, 
and the third due to subtransient effects of the amortisseur. The 
subtransient component of the experimental wave is of very short 
duration as compared with that portion caused by the field transient. 
Likewise, the armature transient is of short duration as compared to 
the field transient. The alternator block diagram as developed in this 
thesis has neglected the effects of the subtransient and armature tran- 
sient components (see Appendix B). A study of the transients due to 
load changes in figures II through X shows the discrepancy between 
computer and experimental curves caused by neglecting these com- 
ponents. Their presence can only be noted for a short interval follow- 
ing the initiation of the transients, and they have completely disappeared 
after 0.1 seconds. After their disappearance, the experimental tran- 
sient closely follows the computer transient. This result justifies 
the assumption that the alternator behaves like a single time constant 
machine, 

In each of the nine transient comparisons the time constant 
associated with the experimental transient was somewhat less than 


the computer time constant. The best explanation for this consistent 
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error is that the measured value of field leakage reactance, Xp: was 
probably in error since this was the least reliable parameter measure- 
ment. It was obtained by first measuring the time constant of the field 
transient associated with a step in field voltage with open circuit ar- 
mature terminals. The time constant was then multiplied by the field 
resistance to obtain the field inductance. The leakage term was ob- 
tained by subtracting the mutual inductance term, common to both 
field and armature, from the value of field inductance. Thus, the 
field leakage inductance was the difference between two values close 
in magnitude and, hence, in subject to considerable error. A smaller 
value of field leakage inductance would decrease the computer time 


constant and make its transient approach the experimental transient. 


The discrepancies occurring between computer and experimental 
transients for field steps were greater than for load changes. This 
was probably due to the method used for applying the field steps. In 
order to avoid changing the field resistance, the step was applied by 
shifting the excitation voltage source of 112 volts D.C. to another 
voltage source of approximately 70 volts D.C. A double throw switch 
was employed for this purpose. The switching delay was about 0.02 
seconds during which time the excitation voltage was zero. Thus, 
during the delay period, the armature voltage was dropping toward 
zero as a final asymptote. The resulting experimental transient thus 
dropped more rapidly than would be the case if no switching delay 


were introduced. 


The block diagram developed in this thesis does not account 


for changes in speed. As was previously explained, the initial or 
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final values of the experimental transients were adjusted in each case 
where they were graphically compared with the computer transients. 
In general, any change in load or excitation will cause a generator 
speed variation. These speed changes should be incorporated into 
the block diagram representation in order to extend its employment 
in the study of the voltage regulation problem. A method to properly 
account for speed in the alternator representation should be a subject 


for further study. 


The final representation was developed assuming balanced load 
conditions. Any study of the voltage regulation problem should con- 
sider unbalanced load conditions due to faults in a single phase of the 
generator load. Thus, an additional recommendation for further in- 
vestigation would be to study the possibility of including unbalanced 


loads in the generator-load representation. 
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V 


CONCLUSIONS 


(1) The block diagram development obtained in this thesis is a valid 
representation of the synchronous alternator. This conclusion is 
based on the close correlation obtained between the computer and 


experimental results in both the transient and steady states. 


(2) The block diagram representation is a comparatively simplified 
form. This simplification results in some slight departures 
from actual steady state values in the saturated regions of the 


magnetization curve. 


(3) At unity power factor, the computer steady state characteristic 
lies above the experimental curve. At zero power factor, the 
computer characteristic lies below the experimental curve. At 
load power factors intermediate between zero and one, the dis- 
crepancy between actual and analytic results would probably be 
less than the maximum error obtained at either 0 P.F. or 1.0 P.F. 
Since the normal load ona generator is rarely at either of these 
extremes, the resulting correlation of computer and experimental 
steady state data utilizing the alternator representation herein 
presented should be better than that obtained for either unity or 


zero power factors. 
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(2) 


(3) 


(4) 


VI 
RECOMMENDATIONS 


Investigation of methods to reduce the steady state errors at high 
values of field excitation with both zero and unity power factor 


loads. 
Investigation of the plausibility of including the effects of armature 
transients and amortisseur subtransients in the alternator block 


diagram. 


Study of a method to properly account for speed variations in the 


representation. 


Study of inclusion of unbalanced loads in the generator-load repre- 


sentation. 
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APPENDIX A 


DEVELOPMENT OF THE LINEAR VOLT-AMPERE EQUATIONS 
DESCRIBING THE SYNCHRONOUS ALTERNATOR 


Introduction 


In the development of the volt-ampere equations for the 
synchronous alternator the methods outlined in reference 10 have 
been followed. The equations are first developed assuming motor 


action and utilizing the machine structure of figure I. 


FIGURE (I) 
Physical Structure of Machine 


DAL, 





Figure (I) places the salient field pole structure on the stator 
although in practice the d. c. field is usually placed on the rotor. Since 
relative motion is the only consideration, this representation is valid. 
However, in order to conform to standard synchronous generator prac- 
tice, the equations will later be changed to show generator action and 


the salient field structure on the rotor. 
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The development below also makes the following assumptions: 


(1) Geometrical symmetry exists on at least one side of 
the air gap. 


(2) No slot openings exist. 

(3) No saturation occurs. 

(4) Sinusoidal mutual air gap flux density exists. 

(5) Amortisseur windings are neglected. 
Assumption (1) is valid since the machine tested has a cylindrical stator 
structure. Assumptions (2) and (3) arc approximations which will be 
used in this portion of the report in order to develop linear differential 
equations with constant coefficients. Appendix B describes the method 
for the introduction of the nonlinear effects of saturation. Assumption 
(4) is a series approximation which is used even though the rotor of 
the machine tested is a salient member; this approximation can only 
be justified by the results obtained. For the machine tested, assump- 
tion (5) is not valid and is only used in order to simplify the analysis. 
In order to include these windings in the analysis it is necessary to add 


a second winding on the stator member of figure (I). 


The development to be presented also assumes that the machine 
equations can be derived from linear lumped circuits in relative motion. 
Therefore, the three phase belts of the armature are replaced by three 
coils each of N equivalent turns which are sinusoidally distributed in 
space. Matrix algebra is used to aid in the development of the equa- 


tions and reference 10 has been used throughout as a guide. 


Phase Components 


The electrical equations of motion for the salient pole machine 
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of figure (I) in terms of the phase variables are 
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Finally, in terms of the phase components the equations become 


mr 
VN R +p bath, cos 29] p FL, +1, cos(20 + d)} | 
r 
Vy |- lp 上 Lu+L。cos(20+d| — R"+p|[L. +L. eos(20 + 20) | 
E Eius Been] o8 | 
T M cos f M cos( +d) | 
' — 
'P [-L +L, cos(20 + 2d) M cos Ø I, 
ED ÉL, +L, cos 20] M cos(@ + d) I, (13) 
r 
centen) Messa) |L, 
M cos(@ + 2d) Re +p L; I, 


Symmetrical Components 
In order to simplfy the equations, the Lyon transformation is used 
to obtain symmetrical components. To accomplish this simplification, 


use is made of the transformations 


1 1 1 
Qu - = loa = where a = e312 713) (14) 
l 2 a 


ei 





and 


0 = 1 ae a (15) 


and the voltages, currents, and impedance are converted to symmetrical 


components as follows: 
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Combining equations (16) through (26) in order to obtain a complete 


volt-ampere equation in Lyon instantaneous symmetrical components, 


A R'+p(L. -2L,) 0 0 á I 
Xl 0 R'+p(L +L.) pele 3% | pM eid || (27 
ME nh» | B zh, o? a Í R'+ p(L,+ Ly)' ‘> el? | [1 
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Rotating Components 


In order to obtain linear differential equations with constant 
coefficients the positional dependency must be removed and a speed 
dependency substituted. Then, if constant speed is assumed, the 
equations become linear. For a salient pole machine the necessary 
transformations can only be made by referring the symmetrical mem- 
ber to the salient member; in this instance this will mean that the rotor 
quantities will be referred to the stator. In equation (27) only the +, 

-, and field components are positionally dependent and torque-producing; 
therefore, only these components need to be transformed. In order to 
convert symmetrical components to rotating components the transfor - 


mations shown in equation (28) will be used. 


eJ o 9 EJE A 
B=|o ed? o; @'-lo ed? o w) 
0 O l 0 A 


With these transformations the voltages, currents, and impendances 


in rotating components become 
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After completing the necessary matrix multiplications the completed 


equations in rotating components become 
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Park's Direct and Quadrature Two-Axis Components 

Since the rotating components are speed dependent with rotor 
quantities referred to the stator, the real direct and quadrature two- 
axis components can be obtained directly from the complex rotating 


components by using the change of variable of equation (32). 
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The volt-ampere equation in two axis components becomes 
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Combining the results of equations (34), (35), and (36), the completed 


equations in Park's two-axis components are 
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Defining the direct axis and quadrature axis inductances as 
1 d 
d a b 2 c = 
5 _3 (38) 
2 Lt L, 5 L 
The equations for motor action and the physical structure of 
figure (I) become 
r 
Mo R + p(L, - 2L.) 
Vi = 0 (39) 
V 0 
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Generator Equations 


The volt-ampere equations have been derived for motor action 


assuming clockwise rotation of the armature as indicated in figure (I). 
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In order to obtain generator equations the following changes must be 
made: 

(1) Since the field (rotor) rotates in a clockwise 
direction, the armature must revolve counterclock- 
wise or opposite to the direction shown in figure (1). 
Therefore, all speed terms in the matrix must have 
their signs changed. 

(2) The matrix as given is developed in terms of voltage 
drops. In order for the equations to be in terms of 
generator voltage rises, all V terms must have their 
Signs changed, i.e., e = -V. 

(3) The electromagnetic torque for a motor is positive 
while for a generator it is negative. Since the sign 
of the torque is governed by the polarity of the exci- 
tation current, I. must be made negative for a gen- 
erator. 

Making the changes indicated above, the volt-ampere matrix. for a 


generator becomes 





r 
ESSE DIE. -2L 0 E EM s 
p( a b) | 
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0 -ÖL -R'-pLa Me 
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Thus, the linear volt-ampere equations describing the synchronous 


alternator in Park's d-q two-axis components become 


一 、 r 4 
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APPENDIX B 


DEVELOPMENT OF A COMPUTER BLOCK DIAGRAM 
REPRESENTATION OF THE SYNCHRONOUS ALTERNATOR AND 
ITS CONNECTED LOAD, INCLUDING EFFECTS OF NONLINEARITIES 


Development of Linear Block Diagram 
The linear differential equations derived in Appendix A for a 
salient pole synchronous generator with resistive-inductive loading 
are 
eq = -(R+pL4) +0L Iç t M I; (1) 


eg = OLgly-(R+ pL + MOI, (2) 





FIGURE (1) | 
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These equations can be represented by a block diagram in which er is 
the exciting signal and the direct and quadrature currents and voltages 


feed a load as yet undefined. 


The block diagram obtained above is derived from equations 
which do not GO for the subtransient effects due to the amortis- 
seur. These subtransients are of very short duration compared to 
others which are present, and it is often justified to neglect them in 


considering the problems of voltage regulation and system design. 


It should be noted that ey Eq’ ig, and are actually direct 


q 
current modulating envelopes of an A.C. carrier. Variation in these 


quantities occur only transientwise. 


The equations representing the load to be applied to the generator 
terminals can be obtained in a straightforward manner from a study of 
the vector diagram describing the volt-ampere relationship of the load 
in the steady state. The steady state relationship is not an exact rela- 
tionship, since any changes in the load signal will cause transients to 
Occur. However, these transients, which are unidirectional, do not 
produce appreciable effects and hence are neglected in the same manner 


that the armature transients are neglected as will be shown later. 
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on 
il 


power angle between et 
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absolute values of load impedance 


FIGURE (II) 
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id = i. sin(6 + 0) (4d) 


also: 
eq 7 €¿sinó = i ZL sin ô (5d) 
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Substituting equation (6) into equation (4), 


E 
d Z., sinô 


Sin(ô + 0) 
L 


“d 
ae [sind cos 0 


Zi Sa y sin@cos 4 
f 
ey 3] 
= — [cos 9+ sin 90 cot 6 (7d) 
z [ 
e ! 
but cot 6 = —2 (8d) 
e 
d 
Therefore 


E. a 
d Zj 


lie 
q 2, cos Š 


っ 
Zr cos 6 


e 
z— cos 9 + = sind (9d) 
L 


i. cos (6 + 9) (4q) 


e, cos 0 = i, Z1, cos ô (5q) 


€ 


is - Zj COS 6 (6q) 


に 
cos(6 * 0) 


[cos 6 cos 0 
L - sin 6 sin 9 


e 
. BA — 
E > [cos 9 sin 0 tan 6] (7q) 


ed 
tand = e (8q) 
q 
M “d 
cos 9 - =— sin 9 (9q) 
ー Ar, 


Applying equations (9d) and (9q), the load can now be represented 


in block diagram form and added to the block diagram earlier obtained 


for the generator. 














<ーーーーーーーーー Generator 
FIGURE (Ill) | 


Basic Assumptions 


The block diagram represented above does not account for the 
saturation effects normally encountered in a magnetic circuit. To 
account for the saturation of the machine, the nonlinearity is approxi- 
mated by means of a single characteristic which can be placed in such 


a position in the block diagram as to correct for the linear, 
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characteristic block diagram obtained above. 


In addition to neglecting subtransient components, it can also be 
seen that armature transients can be neglected. These transients are 
largely responsible for a second harmonic wave which is not detect- 
able on the experimental transient current wave forms, and also 
n for the direct current transients of the output wave which vary 
from phase to phase depending upon the time in each cycle at which 
the transient signal is introduced. It would thus be quite difficult to 


compare these D.C. transients experimentally and analytically. 


An acceptable procedure is to neglect them since they、do not 
produce appreciable effects upon the average terminal voltage of the 
three phases. This assumption amounts to dropping the ply DHL 
and the p M terms from the block diagram. 

In a salient pole machine, the reactance of direct axis armature 
reaction, Xd saturates to a much greater extent than the corres- 
ponding quadrature axis reactance, Kga’ T is, therefore, a reasonable 
assumption that X. , and thus La does not saturate, but remains 


q 
linear. 


Two final assumptions are that the armature resistance is 
negligible, and that the speed remains constant corresponding to an 


electrical frequency of 60 cps. 


Per Unit System 


In order to simplify the linear block diagram so as to more easily 


handle saturation effects, a per unit system will be used as here 
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defined. The stator quantities based on machine rating are defined as: 


BASE ARMATURE QUANTITIES 


① 
I! 


rated peak phase voltage 








ao 

— - rated peak phase current 
es = oL. o E E 9 las - armature base reactance 

La L 
X. = NET om 
d Lao q Lao 

Em - Bad 

ad Lo 


The per unit system as now defined for the rotor circuit is a 
reciprocal system in which the base field current is the field current 
which will induce in each stator phase a voltage equal to Xd -— (8) 
The per unit impedances of the elements in the rotor circuit (field 
circuit) are obtained by transferring the corresponding physical im- 
pedances to the stator by an effective stator-rotor turn ratio, and 
dividing the transferred values by Xo This makes X nf = X d = Af 
= X nd in the per unit system. Thus, the following per unit field 


quantities can be defined: 


that field current which is read directly on the air-gap 
line of the no-load saturation curve corresponding to 
a per unit armature voltage equal to per unit X d 


S 
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p = per unit time differential operator = d 


Block Diagram Including Nonlinear Effects 


If the block diagram parameters are expressed in the per unit 
reciprocal system as defined in the preceding section, and the assump- 


tions of Section IJ are applied, equations (1), (2), and (3) can be written. 


ca (10) 
Sonar "US * nlla (11) 
ep = (Re + PAL gly - PX pd la (12) 


The effect of saturation is introduced in a single characteristic 
which effectively converts an open circuit voltage to a linearized air 
gap voltage. This introduction of the nonlinearity is accomplished by 
changing from the open circuit field current to a hypothetical air gap 


current as shown below. 


Q T» 
€; 
| | 
C) | | lito 4 3 | 
| j | 
| | | 
| 
| | | 
| | | 
I 
Sea. oc 
— — 
1。 X. 


FIGURE (IV) 
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The complete block diagram which accounts for saturation of the 


magnetizing circuit (assuming Ka linear) is given below in per unit. 


E E 


+ 
Es 


A 











FIGURE (V) 


The load portion of the block diagram can be further reduced if 
the load parameters are considered constant. For the experimental 
verification of the block diagram, Z1, and Ka are constant, and for a 
specific load condition, @ is constant. Thus, eq = t (eq) and ia = fte) 
for constant power factor. 

Thus: 
id = K1(9)e。 ( 13) 
and 


eq * Ko(9)e。 (14) 
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where: 


x sod p " 0 
BR (ei. = : + ラーー (15) 
1 Zi Zj, + sin orx Z1, 
and 
X cos 0 
K,(9) = 7 (16) 


L + sin NK) 


The per phase terminal voltage e, can be written 


t 


[Kore] den eq. 


P JI. to] F (17) 


Thus, the generator and load can be represented as 


ES = Wed 2 d 
md K, (8) +] 





FIGURE (VI) 


Computer Representation 


The block diagram of figure (VI) is not directly solvable by ` 
analogue computer techniques because of the presence of differentiators. 


The representation of a differentiation on an analogue computer creates 
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an unusual amount of ''noise'' in the solution. It is therefore customary 
to transform a diagram involving differentiators into an equivalent 
diagram employing integrators only. The computer representation 
obtained is given in figure (VII). It should be noted that the saturation 


characteristic used is the inverse of that depicted in figure (VI). 
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FIGURE (VII) 
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NRENDIX C 


EXPERIMENTAL DETERMINATION OF MACHINE PARAMETERS 
AND TRANSIENT PHENOMENA 


Machine Parameters 


The various machine parameters were determined from the 
following tests and measurements: slip test, open circuit test, short 
circuit test, and a measurement of the field time constant. An oscillo- 
graph was used to record both the slip test and field time constant 
measurements. As shown by the calculations below, the direct axis 
synchronous reactance, X,, was found from the short circuit and open 
circuit tests (see figure II). The quadrature axis synchronous reactance, 
Xg was then calculated from the value of Xj already found and from the 


fatio Xa Xd as determined from the slip test (see figure 1). 





FIGURE I. SLIP TEST 
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From the short circuit and open circuit tests, 





El 
Xi - I 28 = . 4 = 3.285 ohms. (1) 
Me 8.233 


From the slip test, 


=e 1 0405.6 


= (2) 
Therefore 


x, = 0.662 X, = (0. 662)(3. 285) = 2.175 ohms. (3) 


The direct and quadrature synchronous reactances, Xq 


can be expressed in per unit as the sum of an armature leakage reac- 


and X , 
q 


tance Xo} and a reactance X d corresponding to the direct axis arma- 


ture reaction flux: 


NX (4) 


al 


tu (5) 


In a similar manner the per unit field self reactance Xr can also be 
expressed as the sum of a field leakage reactance Xi and the field 


magnetizing reactance X which in the reciprocal per unit system 


fm’ 


(8) is equal to Zen 


Kee = Bey + Xe = Mp + ma (6) 


In order to obtain an adequate block diagram representation it 
was necessary to determine the alternator field resistance Rf， a proper 
load representation, a method for insertions into the block diagram the 


nonlinear effects of saturation, and values for Xp, X and Zee which 


al’ 
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have already been defined. For the computation of the values above, 
a reciprocal per unit system (8) was used in which base field current 
is so defined that X q the reactance of direct axis armature reaction, 
is the equivalent per unit magnetizing reactance between armature and 
field, X nd Therefore, base field current can be read directly on the 
air gap line of the no-load saturation curve corresponding to a per unit 


armature voltage equal to per unit X d 


Of the parameters needed for the block diagram, Appendix B 
describes the methods used to determine a load representation and the 
nonlinear saturation characteristic. The calculations below determine 
X 


the remaining parameters: R, X and X nd In all cases per 


fl’ “al” 
unit values are found by the reciprocal system described in reference 


8 and are referred to the armature. 


The magnetizing reactance X nd and the armature leakage 


reactance X] are found from the Potier triangle of figure (III). The 


distance ac is proportional to Xv bc is proportional to X. d - X nd’ and 
ab is proportional to X It follows that 
X 
Eus obe 7895 (7) 
Xa ac 
and 
X 
nu MP. 02105 (8) 
Xq ac 
Using the value for Xd already determined, 
X „q 70.7895 Xj 72.59 ohms (9) 
- a 10 
X,] 5 0.2105 X¿ = 0.692 ohms (10) 
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and in per unit values 


md 


al 


where the base 


e and i are 
ao ao 


Kad q 
a = a T 0. 592 per unit (11) 
ao 
= 
= —— = 0.1576 per unit (12) 
Aao 


armature reactance, ee = em = 4 39 ohmswnd 


Í 
ao ao 


nameplate values (peak). 


The field leakage reactance Xp was obtained by finding the field 


time constant from the application of a step change in field excitation 


and a record of the resulting transient on an oscillogram (see figure IV). 





FIGURE (IV) 


SN 


Thus, knowing the field time constant r, and having measured the field 


resistance Re» the field inductance could be calculated. 


x AE o 6.96 henries (1 3) 


The determination of per unit field self-reactance Ker follows. 


Lid 


3 f fo gv ; 
A = € { J = 0.69 per unit (14) 
E" 2 372 an 
where L o X, o! 9 
lo - peak armature rated current 
ito - base field current read on the 
air gap line of the no-load satura- 
tion curve corresponding to a per 
unit armature voltage equal to per 
unit X d: 
= 1.7 amps 
| 
Since Xr = AH + X nd’ 
Xp = 0.69 - 0.592 = 0.098 per unit (15) 


. The only per unit values remaining to be found for a complete 


block diagram representation can be calculated from the following 


formulas: 
X 
d ohms) 
X = ——NIF-x = 0. 749 per unit (16) 
(per unit) a ! 
Xa 
x = (ohms) . p 49s per unit (17) 
Viper unit) Lao 
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R š 
f i q 
R „ 3 (ohms) T =2.63x 10 * R (18) 


per unit) e Xao ao (ohms) 


e 


I volts) ifo 
= ee | = 0.0173 per unit (19) 
ao ao 


= 
(per unit) 


Transient Phenomena 


The dynamic performance of the synchronous alternator was 

tested at three widely varying load conditions: pure inductive load 

(pF = 0), pure resistive load (pF = 1), and a combination inductive- 
resistive load (pF = 0.752). At each load condition step changes were 
applied to the alternator and the ent effects upon armature ter- 
minal voltage, armature current, and field voltage were recorded by 
means of oscillograms. The step changes consisted of going from full 
load to no load, going from no load to full load, and applying a step in 


the field excitation. 


The oscillograms of figure (V) through (XIII) indicate the steady 
state values before and after each step change and the time constants 
associated with the transient phenomena for the three load conditions 


used. 


53 





Experimental 


U 


NW ! om HN 
MEM MH NETT I iil 








Transients - PF = 0 


Figure V 


Full Load to No Load 


e, = 0.988 per unit 
] 
e, = 1.372 per unit 
La 
T 7 0.10 sec 
Figure VI 


No Load to Full Load 


e, = 1.363 per unit 
l 
, = 0.982 per unit 


2 
T= 0. 125 Sec. 


e 


Figure VII 
Field Step 
e, = 1.0 per unit 


e, = 0.705 per unit 


+ = 0.155 Sec. 





Experimental Transients - PF = 1.0 


Figure VIII 
Full Load to No Load 


e. = 1. O per unit 


① 
I! 


1.32 per unit 


T = 0.125 sec. 


Figure IX 
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Figure XI 
Full Load to No Load 


En: 0.991 per unit 
l 


e, = 1.446 per unit 
E 


- 0.11 sec. 


Figure XII 
No Load to Full Load 


e, 7 1. 423 per unit 


Hl 


e, = 0.978 per unit 


ラ 


+r = 0.14 sec. 


Figure XIII 
Field Step 


e, = 1.009 per unit 


加 
e, =0.737 per unit 
ta 


Tr = 0.145 sec. 





APPENDIX D 


COMPUTER DETERMINATION OF TRANSIENT PHENOMENA 


In order to check the validity of the synchronous generator 
representation developed in Appendix B, it is necessary to carry out 
both experimental and analytic tests which can be correlated in the 
transient and steady states. The purpose of this appendix is to obtain 
the analytic results of various changes in excitation and load conditions 


on the generator output. 


One of the most difficult problems in the analysis and synthesis 
of the synchronous generator arises from consideration of saturation 
and other nonlinearities. The method employed in this thesis consists 
of concentrating all the generator nonlinearities at a single point ina 


block diagram which can be represented on an analog computer. 


The principle of operation of the analog computer is based on the 
fact that analogies can be drawn among the various types of physical 
systems when the differential equations describing the systems are 
comparable. Of the various types of analog computers, the one most 
frequently used is the d-c analog computer in which direct voltages are 
used to represent the dependent variables. The simple computing net- 
works in the d-c analog computer are made up of RLC networks and 


d-c amplifiers. 


Although the analog computer can carry out the differentiation 
operation, the use of differentiators is avoided when possible because 


they tend to amplify signal noise and cause unwanted high frequency 
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oscillations. Usually when differentiators are present in a circuit 
they can be replaced with integrators by proper adjustment of the 


analog circuit. 


There are several ways of representing system nonlinearities on 
the analog computer. The one resorted to in the solution of the syn- 
chronous generator problem is to construct a nonlinear characteristic 
curve relating the field current oe) which causes an open circuit 
armature voltage, toa field current, La y) which generates the same 
armature voltage on the air gap line of the open circuit characteristic 
curve. From the curve obtained, a photographic plate is prepared for 
use in the computer function generator. A signal input to the function 
generator causes a cathode ray beam to ride on the photographic con- 
tour generating a continuous output signal. Figure (I) shows the photo- 


graphic plate used in the computer work described in this appendix. 





FIGURE (i). COMPUTER FUNCTION GENERATOR PLATE 
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The block diagram of figure (VII), appendix B, must be 
represented on the analog computer in terms of the computer tools 
available; i.e., integrators, summers, potentiometers, voltages, etc. 


This representation is given in figure (II). 


A shift in the time scale of the integrators was performed by 
dividing all integrator input potentiometer settings by the time scale, 
aj: Thus, if a, were 100, the computer time to perform a given opera- 
tion would be one hundred times as long as for the real time alternator 
operation. This shift in time scale was made necessary for two reasons: 
to allow the transients to be recorded at a sufficiently slow speed that 
the dynamic response of the recording instruments would not be ex- 
ceeded, and also to adjust the signal level of the integrator output to 


insure that the maximum signal level of 100 volts would not be exceeded. 


The differential operator for this per unit work is taken equal to 
d/d(w t) where oj is the electrical frequency of the generator, 12077 
radians per second. Thus, p = (1/1207T) p' where p! is d/dt. The 
integration term, 1/p, therefore equals (12077)(1/p'). In order to 
reduce the per unit block diagram integration to a pure integration, the 
integrator output signal is then divided by 12077. The RC network 
shown at the input No. 8 summer is a stabilization network used to 
prevent system oscillation. 

Transients are initiated by throwing a switch on the computer. 


The switch movements cause a step in the field voltage or a change in 


the load condition between full load and no load. 


A preliminary step in the use of the computer is to check the 


steady-state values of field current vs. armature voltage for various 


Ln 
sO 
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conditions of loading. This step consisted of obtaining the generator 
saturation curves for no load, 1.0 power factor, and zero power fac- 


tor. These curves are given in figure (III). 


The conditions under which the experimental work of appendix C 
was determined were duplicated on the computer. The computer poten- 
tiometer settings which were used for each solution are tabulated in 
table (I). The transient solutions obtained, with their steady-state and 


transient values, are pictured in figures (IV) through (XII). 
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APPENDIX E 


SUMMARY OF DATA AND CALCULATIONS 


Definition of Parameters 


ao 


al 


ps 


ad 
Xaq 


md 


field voltage 


direct axis armature voltage 


quadrature axis armature voltage 


peak armature rated voltage, line to line 


E + ea = armature terminal voltage 


field current 


direct axis armature current 


= quadrature axis armature current 


peak armature rated current, per phase 
"^ 
!d 
field self reactance 


+ io = armature terminal current, per phase 


= field magnetizing reactance 
= field leakage reactance 


= direct axis synchronous reactance 


quadrature axis synchronous reactance 


= I = base armature reactance 
2 ao INC 


armature leakage reactance 


= reactance of the direct axis armature reaction 


reactance of the quadrature axis armature reaction 


magnetizing reactance between armature and field 


mutual inductance between armature and field 


direct axis synchronous inductance 
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Le = quadrature axis synchronous inductance 
Re = field resistance 
Tr = open circuit field time constant 
Zt, = load impedance 
wo = rated electrical frequency 
Summary 
Parameter Actual values : Per unit values 
e 220 |2 volts 
ao 
T 28.9 (2 amps 
X 4.39 ohms 
ao 
0.01164 henries 
ao 
Wo 120 て rad/sec 1.0 
6+ 112 volts 0.0173 
X1 3.285 ohms 0.749 
Xa 2.175 ohms 0.495 
X d 2.59 ohms 0:592 
X nd 050592 
X 0.692 ohms 0.1576 
Xe 0.098 
T+ 0.0417 sec 
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APPENDIX F 


EQUIPMENT 


Generator 
Westinghouse synchronous alternator, salient pole (rotor), 


with amortisseur windings. 


KVA 11 h pF 90 

Volts 220 Excitation Volts 125 
Amps 28.9 Excitation Amps 5 

Cycles 60 Style No. 897363 
RPM 1200 Serial No. 1/ 5c3982 


Driving Motor 


Westinghouse direct current motor, shunt wound. 


Volts 230 serial No. 955771 

Open Intermittent H.P. 15 Type RF 10 

Speed 550/1650 Form A17 
Computer 


Reeves Electronic Analog Computer (REAC). 
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